Delocalized quasiparticles throughout the vortex state in s-wave superconductor 
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Quasiparticle transport in the vortex state of an s-wave superconductor at T —* was investi- 
gated by measuring the thermal conductivity of LUM2B2C down to 70 mK in a magnetic field 
perpendicular to the heat current. In zero field, there is no electronic conduction, as expected for a 
superconducting gap without nodes. However, as soon as vortices enter the sample quasiparticles are 
seen to conduct remarkably well, even better than they would in a typical d-wave superconductor. 
This is in stark conflict with the widely held view that quasiparticle states in s-wave superconductors 
just above H c \ should be localized and bound to the vortex core. 

PACS numbers: 74.70.Tx, 74.25.Fy 



The physics of weakly interacting electrons in non- 
magnetic weakly disordered metals and superconductors 
with three-dimensional crystal structures is thought to 
be extremely well understood, in the following limits: 
the s-wave superconductor in zero magnetic field and 
the metallic state at low temperature. The former is 
described by the theory of Bardeen, Cooper and Schri- 
effer (BCS) and the latter by the Fermi-liquid theory of 
Landau, two cornerstones of solid state theory. These 
two limits can be achieved in the same material and the 
electron system can be made to go smoothly from one to 
the other simply by applying a magnetic field H , going 
from zero to above H C 2, the upper critical field beyond 
which superconductivity is destroyed. 

The s-wave superconductor at T — > and H = is 
a thermal insulator, in that the electron system cannot 
transport any heat: the Cooper pairs carry no entropy 
and there are no quasiparticles. On the other hand, the 
normal state above H C 2 is a thermal metal. The two 
states may be characterized by their heat conductivity 
k(T), through the value of the residual linear term kq/T, 
defined as the limit of k/T when T — > 0. kq/T is zero in 
the s-wave superconductor and finite in the metal, where 
it is given precisely by the Wiedemann-Franz law. This 
law is a universal relation between charge conductivity 
and heat conductivity at T — > 0, given by kq/T = Lq/ po, 
where pa is the residual electrical resistivity. 

A basic question arises: how does the electron sys- 
tem evolve between these two well-defined and well- 
understood limits, these two distinct states of matter? 
Remarkably, this question is still largely unanswered. 

In superconductors of the second kind (type II), the 
electron system responds to an applied magnetic field 
greater than a critical value H c \ by generating vortices, 
i.e. lines of quantized flux at the core of which the 



superconducting gap goes to zero. Our basic question 
then becomes: what is the nature of the electronic states 
as a function of vortex density? The density is given 
strictly by the magnetic field strength B as &o/B, where 
$0 = hc/2e is the quantum of flux. Considerable progress 
has been made over the past decade in elucidating the ex- 
treme limits of very low density and very high density. A 
single isolated vortex in an s-wave superconductor with 
an isotropic zero-field gap Ao introduces a discrete set of 
energy levels inside the vortex core, the lowest level lying 
an energy Ao/2fci?£ above the Fermi energy, with higher 
levels spaced equally by twice that amount pj, where 
kp is the Fermi wavevector and £ is the vortex core ra- 
dius. The density of states associated with these levels 
is roughly that of the normal state, but limited to an 
area 7r£ 2 , where £ is of order 100-1000 A in most s-wave 
superconductors. The important point is that the states 
associated with a single vortex are not extended but lo- 
calized, bound to the vortex core. The detailed variation 
of the local density of states associated with a vortex was 
mapped out using STM on NbSe2 by Hess et al. || and 
interpreted in terms of states bound to the vortex core 
||. In the presence of a vortex lattice, the quasiparticle 
states form bands ^||| , but at fields low enough that the 
intervortex separation d is much larger than the vortex 
core size, the low- lying bands are very flat |j| , essentially 
identical to the localized levels of the isolated vortex (and 
therefore also gapped). In other words, when the ratio 
£/<i is very small, the tunneling between states on neigh- 
bouring vortices is exponentially small, causing very little 
dispersion. 

In the other limit of fields close to H C 2 , quasiparticles 
are delocalized, as inferred for example from the obser- 
vation of de Haas-van Alphen oscillations in type-II su- 
perconductors such as NbSe 2 @, V 3 Si and YNi 2 B 2 C 
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[gj. In this limit, calculated bands are broadened Landau 
levels with zero-energy quasiparticles occuring at several 
k points in the magnetic Brillouin zone [Q. Tesanovic 
and Sacramento have argued that the high-held gapless 
phase is separated from the low-field gapped phase by a 
quantum level crossing transition Q . 

The picture for very low and very high vortex densities 
is therefore expected to be much the same as at H = 
and H = H C 2, namely that of localized vs delocalized 
quasiparticles, respectively. The best way to investigate 
the nature of quasiparticle states and the possibility of a 
metal-insulator transition is to measure heat transport. 
Note that charge transport is much less ideal, either the- 
oretically since charge is not a conserved quantity, or ex- 
perimentally since charge dynamics is dominated by the 
dissipation caused by moving vortices (flux flow regime) . 
In this Letter, we report on the first systematic study of 
heat transport throughout the vortex state of an s-wave 
superconductor. Contrary to all expectations, we observe 
highly delocalized quasiparticles at temperatures as low 
as T c /200 for all values of the field in the vortex state. 
There is no indication of an insulator-to-metal transition 
between H c \ and H C 2, eventhough the ratio <i/£ starts off 
as large as 10-20 at H c \. 

The borocarbide material LUNI2B2C was chosen be- 
cause of its large critical field ratio H^jHd — 100, which 
allows access to a wide range of intervortex separations. 
It is a non-magnetic superconductor thought to be char- 
acterized by s-wave pairing and a gap without nodes. 
Its crystal structure is tetragonal and its Fermi surface 
properties nearly isotropic. The thermal conductivity 
was measured in a dilution refrigerator using a standard 
steady-state technique. The temperature was swept at 
fixed magnetic fields, from 70 mK up, in fields ranging 
from to 8 T. The field was applied perpendicular to the 
heat current and parallel to the c-axis ([001] direction). 
Below 1.5 T, the sample was cooled in the field to ensure 
a good field homogeneity and a well-ordered vortex lat- 
tice. Above 1.5 T, the results were independent of cooling 
procedure. The single crystal was grown by a melting 
flux method described elsewhere |10| . The sample was 
shaped and polished into a rectangular parallelepiped of 
width 0.495 mm (along [010]) and thickness 0.233 mm 
(along [001]), with a 1.59 mm separation between con- 
tacts (along [100]). The same contacts were used for 
electrical resistivity and thermal conductivity measure- 
ments, made with non-superconducting solder and yield- 
ing contact resistances lower than 5 mil at 300 K. In zero 
magnetic field, the resistivity was /?(300 K) = 35 /ifi cm 
at room temperature and pa ~ 1.30 pfl cm at T — ► 0, 
using an extrapolation of the form p = po + AT 2 between 
T c and 50 K, yielding an RRR of 27. T c was found to 
be 16 K with a transition width of 1 K. Magnetoresistive 
measurements performed at 1.8 K yielded H C 2 = 7 T and 
p = 1.67 pfl cm just above H C 2- The zero-temperature 
coherence length £ = 70 A, using H c2 (0) 
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FIG. 1. Thermal conductivity of LuNi2B2C between 2.5 
and 25 K. The superconducting transition occurs at T c = 
16 K and can be seen as a kink in the thermal conductivity. 
The large peak centered at 5 K is due to a rapid increase in 
the phonon mean free path as the electrons condense below 
T c . Inset: Thermal conductivity at T = 2.1 K as a function 
of magnetic field. The downward deviation from a constant 
conductivity at a field H c i ~ 60 mT signals the onset of vortex 
entry into the sample. 

Note that the radius of the vortex cores are somewhat 
larger than this at low fields (by a factor 2 or so) as mea- 
sured by muon spin rotation (//SR) in this compound pl| . 
A similar expansion was previously reported in NbSc2 
& 

The thermal conductivity of LuNi2B2C in zero mag- 
netic field is plotted in Fig. 1 as k(T) vs T up to T c 
and beyond. The total conductivity is the sum of an 
electronic and a phononic contribution: k = n e + n p h. 
Above T c , k is dominated by electrons, as n p h is severely 
limited by electron scattering. As the temperature is 
decreased below T c , two things happen as electrons con- 
dense into Cooper pairs: K e decreases as there are fewer 
and fewer thermally excited quasiparticles and K p h in- 
creases as fewer unpaired electrons are available to scat- 
ter phonons. The first effect is seen as a dip right at 
T c and the latter leads to a huge increase in phonon 
mean free path giving rise to the pronounced peak at 
5 K. This increase is expected to stop when the num- 
ber of thermally excited quasiparticles becomes vanish- 
ingly small. In a typical BCS superconductor such as 
Nb, the peak occurs roughly at T c /5, the temperature at 
which n e becomes negligible (less than 1-2% of its nor- 
mal state value). The fact that in LuNi 2 B 2 C the peak 
is even higher in temperature is an indication that there 
are very few quasiparticles below T c /3, which tends to 
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confirm the s-wave nature of the gap and rule out any 
sizable gap anisotropy. At temperatures below T c /3, k 
is essentially all due to phonons, whose mean free path 
eventually grows to reach the size of the crystal, at which 
point a cubic temperature dependence is observed. This 
is shown in Fig. 2, in a plot of k/T vs T 2 , which al- 
lows us to easily separate the electronic term linear in 
T from the phononic term cubic in T. The zero-field 
curve is indeed linear in such a plot (below 150 mK), 
with zero intercept and a slope in quantitative agree- 
ment with the known sound velocities and sample di- 
mensions, as reported earlier |l3| ]. The absence of any 
residual linear term is a second confirmation that there 
are no nodes in the gap (at least in the basal plane). 
To be more specific, the data sets an upper limit of 
k /T ~ 0.01 mW K~ 2 cm" 1 ~ (k n /T) (1/2000) on a pos- 
sible residual linear term arising from nodes in the basal 
plane. Theoretically, the value expected from a line node 
in the basal plane is roughly k /T ~ (kn /T)(KT/2Ao), 
where T = \/2tn is the impurity scattering rate and 
Ao is the gap maximum fl4| . We can estimate rjy, the 
transport relaxation time in the normal state, from ki- 
netic theory: kjv/T = Lo/po = ^-fNVpTN, where 
is the Sommerfeld coefficient of the specific heat and 
vf is the Fermi velocity. The measured values are 
77V = 19.5 mJ FT 2 mole- 1 |J], A = 3.2 k B T c |J| and 
vf — 2 x 10 5 m/s jlT]]. These numbers give a mean free 
path In — vpT ~ 500 A, and ?ir ~ 0.5Ao, from which 
we would expect kq/T ~ jKn/T if there were nodes in 
the gap. This huge discrepancy, by a factor of nearly 
three orders of magnitude, argues strongly against nodes 
in the gap function of LuN^E^C (and the insensitivity 
to impurities argues against an unconventional order pa- 
rameter). 

As seen in Fig. 2, a residual linear term is induced by a 
magnetic field applied perpendicular to the heat current. 
When the field exceeds -ff C 2(0) = 7 T, the conductivity 
is dominated by electrons and k/T ~ K e /T is essentially 
constant, as seen for 8 T, with a value in perfect agree- 
ment with the Wiedemann- Franz law: K e /T = L a /p, 
where p = p{8 T). Turning to very low fields, we use 
the fact that phonons are scattered by vortices to deter- 
mine in situ the field H c \ at which vortices first enter the 
sample. This is shown in the inset of Fig. 1, where n p h 
at T — 2.1 K is seen to suddenly deviate from its field- 
independent value at H c i ~ 60 mT. Note that at 60 mT, 
it is known from small-angle neutron scattering fl8|| that 
vortices form a well-ordered square lattice in LUM2B2C. 
The question of fundamental interest is: how is the elec- 
tronic conduction modified by the presence of a vortex 
lattice? The low field behavior is shown in the lower 
panel of Fig. 2. The curves of k/T vs T 2 for fixed fields 
of 50, 75, 100, 200 and 300 mT are seen to be rigidly 
shifted relative to the zero field curve. The growth of 
the corresponding residual linear term, which is the elec- 
tronic contribution, is plotted in Fig. 3. The striking 
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FIG. 2. Temperature dependence of the thermal conductiv- 
ity at several applied fields in the range 70 mK < T < 300 mK. 
The data is plotted as k/T vs T 2 to separate the electronic 
and phononic contributions. Fig 1 (a) shows the data for H 
= 0, 0.3, 1.5, 4, 6 and 8 T, in increasing order. The solid 
line indicates the value expected from the Wiedemann-Franz 
law above H C 2- Fig. 1 (b) shows the data for H = 0, 50, 75, 
100, 200 and 300 mT, in increasing order. At low fields, the 
phononic term remains largely unchanged as the electronic 
term grows with increasing field. 

feature is that this growth starts right at H c \ , as can be 
seen from the inset. Whereas quasiparticles in an s-wave 
superconductor close to H c i are theoretically expected to 
be localized and thus incapable of transporting heat at 
T — ► 0, the observed quasiparticle conduction is excel- 
lent. LuNi2B2C is therefore a thermal metal (with fully 
delocalized quasiparticles) throughout the vortex state, 
at least down to T c /200. This provides a natural expla- 
nation for the fact that de Haas-van Alphen oscillations 
in YNi2B2C persist down to unusually low fields (H^/b) 
[|| and no localized core states have been detected in 
STM spectroscopy on borocarbide superconductors. 

The growth of quasiparticle conduction with field is 
remarkably fast. At low fields, it is linear in field (see 
inset of Fig. 3) and given by 
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P0 Hc2 

where po is the zero-field normal-state resistivity. In 
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some sense, this suggests that a quasiparticle in the vor- 
tex state conducts heat just as well as in the normal 
state (where ^ = if the density of quasiparticles is 
assumed to grow roughly as B/H C 2- 

It is striking to note that quasiparticles in this s-wave 
superconductor conduct better than they do in a typi- 
cal d-wave superconductor. Indeed, from calculations by 
Vekhter and Houghton [19], 



(2) 



for all H in the vortex state of a d-wave superconductor. 
Note however that the H dependence is different: it is 
sublinear for d-wave at low H, largely due to the \[H 
dependence of the density of states, whereas we find a 
linear dependence at the lowest fields, perhaps even going 
to super linear (see Fig. 3). 

In trying to explain the surprisingly good conduction, 
it is natural to invoke a highly anisotropic gap func- 
tion, and ascribe the excellent heat transport to the field- 
induced excitation of quasiparticles across a very small 
gap in the basal plane. Izawa et al. have recently in- 
terpreted the \fH dependence of the specific heat in 
YM2B2C [0 in terms of a Doppler shift of the quasi- 
particle spectrum akin to what is seen in d-wave super- 
conductors (H]]2l|] but applied in this case to a highly 
anisotropic s-wave gap, with a small minimum gap A m i n . 
This scenario is unsatisfactory in several respects. First, 
it requires a huge gap anisotropy in a Fermi liquid that is 
otherwise rather isotropic. Indeed, because quasiparticle 
conduction starts right at H c \ , the average Doppler shift 
energy E H 



than the minimum gap: 



2lfl at H c i must be greater 



A min < E H (H cl ) ~ A ^/H cl /H c2 = Ao/10 



(3) 



This factor of 10 or more in gap anisotropy should 
be compared to the very small mass tensor anisotropy 
(e.g. from H C 2), equal to 1.01 in YN12B2C and 1.16 in 
LUN12B2C |T^| . Secondly, as mentioned above, it is hard 
to see how phonon conductivity could rise to a huge peak 
at 5 K if there was significant thermal excitation of quasi- 
particles for T < 5 K, as there would be with A m i n ~ 1 K. 
The third problem is that we observe the same excellent 
quasiparticle transport at very low fields in a crystal of 
LUM2B2C doped with 9% Co impurities. The doping 
causes the scattering rate to be 7 times larger than in 
the pure crystal, so that HT ~ 3.5A and In — £• Such 
strong impurity scattering should essentially wipe out all 
anisotropy in the gap function |2^], yet we observe the 
same transport behaviour in the vortex state of both pure 
and impure samples. 

This last fact also argues against the speculation by 
Nohara et al. |23| that the field dependence of the specific 
heat in YM2B2C goes from j(H) ~ y/H to linear in H 
upon doping with 20% Pt because quasiparticles go 




FIG. 3. Magnetic field dependence of the electronic ther- 
mal conductivity k/T at T — > 0. The current is along the 
a-axis [100] of the crystal and perpendicular to the applied 
field. The arrow gives the value of H C 2 obtained from resis- 
tivity. Note the rapid rise in k/T and the saturation above 
H C 2 at a value exactly given by the Wiedemann-Franz law. 
Inset: zoom at low fields, where a linear growth is seen to 
extrapolate precisely at the independently measured H c \. 

from being delocalized (coherent hopping between ad- 
jacent vortices) to being localized (incoherent hopping). 
We are therefore led to regard the presence of delocalized 
quasiparticles as a property of the vortex lattice state and 
not as the result of nodes or deep minima in the gap func- 
tion. To the best of our knowledge, a proper calculation 
of the thermal conductivity in the vortex lattice state of 
an s-wave superconductor at low fields has never been 
reported. It would have to start from the quasiparticle 
band structure of the vortex lattice, such as calculated 
by Yasui and Kita ||. For an isotropic s-wave gap and 
a hexagonal vortex lattice at B = Q.1H C 2, these authors 
find the two lowest energy bands to be at an energy of 
Ai ~ Ao/2fci?£ above Ep and have very low dispersion, 
with a band velocity no greater than v ~ (<zAo/7i7r)/100. 
These two factors conspire to reduce the expected heat 
conduction relative to the normal state very considerably, 
and it may be difficult to reconcile our results with such 
an energy spectrum. First, the existence of a (mini)gap 
Ai should give rise to an exponential decrease of k/T 
with temperature for fc^T < Ai. A rough estimate gives 
Ai ~ 300 mK, using kp ~ 7r/2a and £ ~ 100 A, yet we 
observe no decrease down to 70 mK. 

Even if we assume that the lowest bands are somehow 
thermally populated, it is still difficult to see how they 
could lead to such excellent conduction, given how flat 
they are. Conductivity is weighted by the square of the 
carrier velocity, or band dispersion, which alone should 
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suppress k/T by a factor (v/vf) 2 relative to the normal 
state. This implies a reduction by a factor 2 x 10 7 at 
B = 0.1H C 2- At that field, K e is only down by a factor 
10. It remains to be calculated whether the transport re- 
laxation time t in the vortex state is sufficiently enhanced 
to compensate for that factor of 10 6 . 

Theoretical expectations of a quantum phase transi- 
tion from thermal metal to thermal insulator as the vor- 
tex separation is increased by going to lower fields are 
not confirmed in our experiments. The appearance of a 
gap at a quantum level crossing transition, proposed by 
Tesanovic and Sacramento ||, could conceivably have 
been missed if one postulates that the minigap is in fact 
smaller than our estimate of Ai ~ 300 mK, by a factor 
of 4 (or more). The localization of quasiparticles as a 
result of disorder, discussed by Vishveshwara et al. 
may simply not occur in the samples we studied (both 
pure and Co-doped) because the level of disorder is not 
high enough. 

In summary, highly delocalized quasiparticles are 
present in LUNI2B2C in magnetic fields as low as H c ±, 
in sharp contrast with the standard picture of quasipar- 
ticle states bound to the vortex cores at very low fields. 
We believe that this remarkably good quasiparticle con- 
duction is an intrinsic property of the vortex lattice state 
as opposed to being caused by a hugely reduced gap in 
the basal plane. Calculations of heat transport in the 
vortex lattice state are needed to explore the apparent 
compensation between the very low band dispersion (or 
quasiparticle velocity) and the transport relaxation time. 
Measurements on other superconductors are necessary in 
order to determine whether our findings are peculiar to 
borocarbides or a more generic property of s-wave super- 
conductors. 

We are greatful to Z. Tesanovic, S. Vishveshwara, T. 
Senthil, M.P.A. Fisher and Michael Walker for stim- 
ulating discussions. This work was supported by the 
Canadian Institute for Advanced Research and funded 
by NSERC of Canada. E. B. and C. L. acknowledge 
the support of Walter C. Sumner Fellowships. L. T. ac- 
knowledges the support of a Premier's Research Excel- 
lence Award from the Government of Ontario. Ames lab- 
oratory is operated for the U.S. Department of Energy by 
Iowa State Univesity under contract No. W-7405-ENG- 
82. 



[8] 
[9] 

[10; 

[11 

[12 

[13 

[14 
[15 

[16 
[17 

[18 
[19 

[20 
[21 

[22 

[23 

[24 



Schluter, Phys. Rev. Lett. 65, 1820 (1990). 

M.R. Norman, A.H. MacDonald, and H. Akera, Phys. 

Rev. B 51, 5927 (1995). 

K. Yasui and T. Kita, Phys. Rev. Lett. 83, 4168 (1999). 
J.E. Graebner and M. Robbins, Phys. Rev. Lett. 36, 422 
(1976). 

F. Mueller, D.H. Lowndes, Y.K. Chang, A.J. Arko, and 
R.S. List, Phys. Rev. Lett. 68 3928 (1992); R. Corcoran, 
N. Harrison, S.M. Hayden, P. Meeson, M. Springford, 
and P.J. van der Wei, Phys. Rev. Lett. 72, 701 (1994). 
M. Heinecke and K. Winzer, Z. Phys. B 98, 147 (1995). 
Z. Tesanovic and P.D. Sacramento, Phys. Rev. Lett. 80, 
1521 (1998). 

P.C. Canfield, P.L. Gammel, and D.J. Bishop, Physics 
Today, October 1998, p. 40. 

A. Price, R. Miller et ai, Private Communication. 
J.E. Sonier et al, Phys. Rev. Lett. 79, 1742 (1997). 
E. Boaknin, R.W. Hill, C. Lupien, L. Taillefer, and P.C. 
Canfield, Physica C 341-348, 1845 (2000). 
M. Graf et ai, Phys. Rev. B 53, 15 147 (1996). 

H. Michor, T. Holubar, C. Dusekand, and G. Hilscher, 
Phys. Rev. B 52, 16 165 (1995). 

Y. De Wilde, et al, Phys. Rev. Lett. 78, 4273 (1997). 
V.G. Kogan et al, Phys. Rev. B 55, R8693 (1997); S.B. 
Dugdale et al, Phys. Rev. Lett. 83, 4824 (1999). 
P.L. Gammel et al, Phys. Rev. Lett. 82, 4082 (1999). 

I. Vekhter and A. Houghton, Phys. Rev. Lett. 83, 4626 
(1999). 



K. Izawa et al, cond-mat/0012155 . 

C. Kiibert and P.J. Hirschfeld, Phys. Rev. Lett. 80, 4963 
(1998). 

L.S. Borkowski and P.J. Hirschfeld, Phys. Rev. B 49, 15 
404 (1994). 

M. Nohara, M. Isshiki, F. Sakai, and H. Takagi, cond 



mat/9902264 



S. Vishveshwara, T. Senthil, 
Rev. B 61, 6966 (2000). 



and M.P.A. Fisher, Phys. 



[1] C. Caroli, P.G. de Gennes, and J. Matricon, Phys. Lett. 
9, 307 (1964). 

[2] H.F. Hess, R.B. Robinson, R.C. Dynes, J.M. Valles Jr., 
and J.V. Waszczak, Phys. Rev. Lett. 62, 214 (1989). 

[3] J.D. Shore, M. Huang, AT. Dorsey, and J.P. Sethna, 
Phys. Rev. Lett. 62, 3089 (1989); F. Gygi and M. 



5 



